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Abstract.  This  paper  is  based  on  the  general  thermodynamical  theory  of 
a  Cosserat  continuum  developed  by  Green  and  Laws .  We  present  here 
specific  constitutive  equations  for  a  linear  viscoelastic  material. 

When  the  form  of  the  free  energy  is  restricted  by  certain  symmetry 
conditions ,  the  basic  equations  separate  into  four  groups ,  two  for 
flexure,  one  for  torsion,  and  one  for  extension  of  the  rod.  Thermal 
effects  occur  only  in  the  last  group.  Flexural  and  torsional  wave 
propagation  nlong  an  infinite  rod  are  considered. 


1.  Introduction 

This  paper  presents  a  theory  of  a  one  dimensional,  viscoelastic, 
Cosserat  continuum.  Following  Green  and  Lavs  [ l3 ,  we  shall  define  a  rod 
as  a  one  dimensional  Cosserat  continuum.  Green  and  Laws  have  developed 
an  exact  thermodynamical  theory  of  rod..,  which  is  not  restricted  to  small 
deformations  or  elastic  rods.  Of  course,  a  theory  of  rods  can  also  be 
constructed  by  considering  the  rod  as  a  three  dimensional  body.  The 
equations  governing  the  rod  are  then  obtained  from  the  three  dimensional 
equations  by  introducing  assumptions  or  expansions  based  on  the  "thinness" 
of  the  rod.  Some  results  which  show  that  the  Cosserat  theory  is  a  natural 
first  approximation  to  the  three  dimensional  problem  have  recently  been 
given  by  Green,  Laws  and  Naghdi  [3].  In  Section  5  of  this  paper,  we  shall 
briefly  discuss  how  solutions  from  the  Cosserat  theory  of  rods  can  be 
compared  with  corresponding  exact  solutions  from  the  three  dimensional 
theory  of  linear  viscoelasticity. 


1. 


Recently  Green,  Laws  and  Naghdi  [2]  have  used  the  basic  theory  of 
[l]  to  derive  a  linear  theory  of  straight  elastic  rods,  although  as  the 
authors  note  their  work  could  be  readily  extended  to  the  case  where  the 
rod  has  an  initial  curvature . 

In  this  paper,  using  a  procedure  similar  to  that  of  [2],  we  derive 
a  linear  theory  of  viscoelastic  rods  which  are  initially  straight.  The 
basic  field  equations  of  the  theory,  which  are  the  same  as  those  given 
in  [2],  are  presented  in  Section  2.  In  Section  3  appropriate  constitu¬ 
tive  equations  for  viscoelastic  rods  are  discussed.  The  discussion  of 
constitutive  equations  and  the  thermodynamics  of  viscoelastic  rods  given 
here  is  analogous  to  the  discussion  of  the  general  linear  viscoelastic 
solid  by  Christensen  and  Naghdi  [4].  We  then  restrict  our  attention  to 
the  case  where  the  rod  possesses  certain  symmetries.  The  resultant 
equations  seem  to  correspond  to  those  of  a  rod  considered  as  an  isotropic, 
viscoelastic,  three  dimensional  member  with  a  cross  section  which  is 
symmetric  about  its  principal  axes.  Under  this  symmetry  restriction  the 
equations  of  the  theory  separate  into  four  groups:  Two  governing  flexure, 
one  governing  torsion,  and  one  governing  longitudinal  motions.  We  note 
that  the  temperature  is  present  in  only  the  last  group. 

In  Section  5  we  use  the  equations  of  Sections  3  and  4  to  consider 
the  propagation  of  flexural  waves  in  an  infinite  rod.  The  solution  for 
an  arbitrary  frequency  is  very  complicated,  and  explicit  results  are 
presented  only  for  the  asymptotic  wave  speeds ,  attenuation  and  ampli¬ 
tudes  as  the  frequency  becomes  very  large.  Hiese  results  generalize 
well  known  results  of  Hunter  [5].  We  also  compare  the  solution  of  the 
problem  of  quasi-static  pure  flexure  using  the  Cosserat  theory  with  the 
solution  of  the  corresponding  problem  using  the  exact  solution  obtained 
from  the  three  dimensional  theory  of  linear  viscoelasticity. 
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Finally  in  Section  6  we  consider  the  torsion  equations.  We  intro¬ 
duce  an  additional  symmetry  restriction  on  the  form  of  the  constitutive 
equations.  Under  this  restriction  our  Cosserat  rod  now  seems  to  correspond 
to  a  right  circular  cylinder.  We  then  study  the  propagation  of  torsional 
waves  through  the  Cosserat  rod  and  compare  our  results  with  those  of 
Berry  [6],  who  solved  the  corresponding  three  dimensional  problem. 
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2.  Linear  Theory  of  Rods 

In  the  basic  theory  of  rods  as  developed  by  Green  and  Laws  Cl], 
a  rod  is  defined  to  be  a  curve  (embedded  in  a  Euclidean  3-space)  at  each 
point  of  which  there  are  two  assigned  directors.  Hie  motion  of  the  rod 
at  time  t  is  then  described  by  the  equations 


r 


r(z,t) 


a 


a,(s.t) 


(2.1) 


where  r  is  the  position  vector  of  the  curve,  a^  (o-l,2)  are  the 
directors  and  z  is  a  convected  coordinate  which  defines  points  on  the 
curve.  We  also  define  the  base  vector  a~  along  the  curve  by 


(2.2) 


and  impose  the  restriction 


*-~l»  £3-*  >  0  •  (2*3) 

In  terms  of  a  and  a_  ,  the  basic  kinematical  quantities  may  be  talten 
as 


(2.4)+ 


+ Latin  indices  take  the  values  1,  2  and  3»  and  Greek  indices  the  values 
1  and  2.  Also,  repeated  indices  imply  the  usual  summation  convention. 


‘id 


a.  •  a 
~i  ~d 


da 

dz 
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The  initial  values  of  the  position  vector  and  the  directors  are 
denoted  by  R  ,  A.  ,  and  A0  .  We  assume  that  the  rod  is  initially 
straight  so  that 

S'  2  *3  ,  ■  6i  j  >  (2-5) 

where  A.  are  independent  of  z  and  6..  is  the  Kronecker  delta.  We 
~i  lj 

restrict  our  attention  to  the  case  when  the  subsequent  displacements  of 
the  rod  are  "small".  To  be  consistent  we  will  also  assume  that  the 
changes  in  the  thermodynamic  variables  such  as  temperature  are  also 
"small".  More  precisely,  we  assume  that 

r-zA3  +  «u  ,  a^-A^cb.  ,  T'  ■  T0  ♦  c  T  ,  (2.6) 

where  T7  is  the  temperature  of  the  body  with  value  Tq  in  the  initial 
undeformed  state,  and  e  is  a  small  non-dimensional  parameter.  All 
forced  and  assigned  loads  are  also  assumed  to  be  of  0(e)  .  By  neglecting 
all  terms  of  0(e  )  and  higher  in  the  field  equations  we  obtain  the 
linearized  field  equations.  After  obtaining  the  approximate  equations,  we 
may  formally  set  e=l  in  the  equations  without  loss  of  generality  with  the 
understanding  that  when  the  displacements,  temperature  changes,  forces,  and 
assigned  forces  are  expressed  in  suitable  non-dimensional  forms  they  are  to 
be  considered  small,  i.e.,  «  1.  Since  this  procedure  is  straightforward, 
only  the  final  results  are  presented  here.  Before  writing  down  the 
basic  equations,  however,  we  note  that  as  a  consequence  of  the  lineari¬ 
zation  all  components  of  vectors  such  as  r  and  a.  are  now  referred  to  the 
fixed  orthonormal  base  vectors  A^,  and  hence,  there  is  no  need  to  distinguish 
between  covariant  and  contravariant  components. 


5. 


It  is  convenient  to  define  measures  of  deformation  by 


or  in  the  linearized  theory 


(2.7) 


V-  .  =  b.  .  +  b  . 

ij  id  di 


Also,  by  (2.4)2  and  (2.6).^,  we  have 


9bi; 


az 


9ui 

3i  az 


(2.8) 


where  are  the  components  of  the  displacement  vector 

u  =  u.  A.  . 

~  1 

The  linearized  basic  field  equations  of  the  theory  are  given  in 
[  2  ] .  Here  we  quote  freely  from  [  2  ]  and  also  record  additional  re¬ 
sults  appropriate  to  the  linear  theory  which  will  be  useful 
later.  The  linearized  equations  of  motion  are 


an.  a  u. 

r-i+pf  =  p  — gi 

1  at2 


(2.9) 


"ag  =  ngcy 


p3~  “g 


(2.10) 


where  p  is  the  initial  mass  per  unit  length  of  the  rod,  n^  are  force 
components,  and  f^  are  the  components  of  the  assigned  force  per  unit 
mrss  on  the  rod.  The  quantities  tt**1  are  defined  by 
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(a .li) 


TTofi~  ^cri.  +  dz  5 

where  are  the  components  of  the  director  forces,  and  q^  are  the 

sum  of  the  assigned  director  force  on  the  rod  and  the  director  inertia 
forces . 

Equations  (2.10)  and  (2.11)  can  be  written  in  a  more  familiar  form 
if  we  introduce  vectors  m  and  g  defined  by 


m=  A 
~  My 


X  q 
~ot 


with  components 


(2.12) 


ml“  p23  *  m2  "  ’  p13  *  m3  "  P12  "  P21 
gl“  **23  ’  ®2  "  q13  ’  g3  "  q12  "  P21 


(2.13) 


It  is  then  easily  seen  that  (2.10)  and  (2.11)  are  equivalent  to 


dm, 

X 

dz 

°2  + 

Pgl 

dnig 

dz~  + 

nl  ' 

pg2 

dnu 
dz  + 

Pg3 

-  0 

(2.14) 


and 
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•or* 


8pn 

"ll  *  Wll  +  —  • 

9Ppo 

n22  "  pq22  +  “az"  * 


2tt12  "  2tt21  “  P^12+q21^  +  Sz  ^P12  +  P21^ 


We  must  also  introduce  an  explicit  form  for  the  director  inertia 
terms .  If  we  assume  that  the  contribution  of  the  directors  to  the 
kinetic  energy  is  of  the  form 


1 

2 


“l2l 


(2.16) 


where  the  coefficients  and  are  independent  of  t  and  a  super¬ 

posed  dot  denotes  partial  differentiation  with  respect  to  time,  it  fol¬ 
lows  that 


B2b 


a, 


ti 


p  at2 


(no  sum  on  p)  , 


(2.17) 


where  4  .  are  the  assigned  director  forces . 

P1 

The  energy  equation  appropriate  to  the  linear  theory  is 


■pA  -  p ( T ' S  +  TS  )  +  pr  +  |  tt^  \  nB^B3+  +  |  n3  Y 


’03  Y3P‘ 


+  pai  Kai 


3  33 


££  =  o 

dz 


(2.18) 


where 


y  .  !lii 

Yid  at 


(2.19) 
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A  is  the  Helmholtz  free  energy  per  unit  mass,  r  is  the  heat  supply 
per  unit  mass  per  unit  time  and  h  is  the  flux  of  heat  along  the  rod. 
To  complete  the  basic  equations  we  postulate  an  entropy  production  in¬ 
equality  of  the  form 


pT'S  -  pr  + 


&  _  J L2I  *  o 

Bz  T  Jz  " 


If  we  use  (2.18)  to  eliminate  pr  from  (2.20),  we  obtain  the  reduced 
entropy  inequality 


-p(A+  TS  )♦£  Ygg +  \  "3  V33 


*  fai  ^  i 5  0  • 


9. 


(2.20) 


(2.21) 


3  •  Constitutive  Equations  for  a  Viscoelastic  Rod 


To  obtain  a  complete  theory  we  must  specify  appropriate  constitutive 
equations  for  the  free  energy  A  ,  the  force  vector  m  ,  the  director 
forces  p^  ,  the  quantities  tt^  ,  the  entropy  S  ,  and  the  heat  flux  h  . 
Since  we  wish  to  consider  the  case  where  n.  ,  p  .  ,  etc.  are  linear  func- 
tionals  of  the  kinematical  and  thermal  histories,  it  is  sufficient  to 
consider  a  constitutive  equation  for  the  free  energy  of  the  form 


pA 


pt  «t 

*  pAQ  +  J  Dij(t“T)vij('OdT  +  J  \(t-T )t(t )dT 


»t  #»t 


»t  (»t 


rr  $ij(t-T,t-Tl)YiJ(T)T(T)dT  dT| 


i»t  «t 


in  m(t-T,t-n)r(T)T(Tl)dT  dT|  +  J  Ha.(t-T)Kai  dT 


i.t  -t 


III  H«iSj(t-,'t-’1)‘I,l(T>'6j(’)dTd1 


»t  «t 


J  I  dTi 


-00  -00 


»t  «t 


rr  rai(t-T5t-Tl)Kai(T)T(Tl)dT  dTl  , 


(3.1) 


where  as  we  recall  a  superposed  dot  denotes  partial  differentiation  with 
respect  to  the  time  variable  t  .  Without  loss  of  generality,  we  may 
assume  that 


10. 


J 


It  is  not  necessary  to  record  here  specific  constitutive  equations  for 

n.  ,  P  .  ,  tt  „  and  S  ;  it  suffices  to  assume  that  they  are  linear 
1  ’  *01  o$ 

m  •  * 

integral  operators  over  the  histories  of  and  T  ,  but  that 

they  are  independent  of  the  temperature  gradient . 

Substituting  from  (3-1)  into  the  reduced  entropy  inequality  (2 .20), 
we  obtain 


11. 


pt  *  l*k 

[-pS  -  \(0)  -  J  m(0,t-T)T(T)dT  -  J  . (t-T,0)v. . (t )dT 

'  2  j  GcW°>t-T>\iWdT  '  2  j  »aB<°.t-T)i(T)dT 


-00  -  - 

■  2 1  Fviae(t'T’0)K^(T)dT]Y~“(t)  +  [2n' 


<*i 


rc*e' 


pt  .  pt 

■  ^  J  G83kA^0,t"T^Ykje,^T^dT  “  ^  J  ^6q(°j^'“T’)t(t )dT 

'  4  J„  F(»63(t"T’0)i;ai(T)dT  •  top3(°>1  Ve3(t) 

Pt 

+  [n3“  2D33(0)  ’  2  J  G33k£(0jt"T)Yk£(T)dT 

—00 

-  2  J_  *33(0,t-T)T(T)dT  -  2  J4  F0,133(t-T,0)Kal(T)dT]  V33(t) 


Ofl 

,t 


tp«i-  Hdi(0)  ■  J.  Hoiej(0,t'T)l'ej(T)dT  -  J_  F««.jk<0>t-T>vJk(T)dT 
J 


roi(0>t-T)i'(T)dT]Kai(t)  -A-J  (t-T^  (T)dT 


-  f  ^  X(t-T)i(T)dT  -  f  £  Ho.(t-,)ia.(T)dT  .fg(t)lO  , 


where 


(3.4) 


I 


12. 


I»t  <»t 


kl  J  n  aijki<t-T>t-1'>Yij(T)Vi(1')dT  d1> 

■©  •© 


.t  .t 


J  J  h  *ij(t-T't-1l)vij(T)T(Tl)dT 


dll 


-00  -00 


*t  »t 


+  ij  J  it  m(t-T»t-TOT(T)T(Tl)dT  dll 


.t  ,t 


+  IJ  J  £  Wt'T’t‘1)“*i(T)V1')dT  t,1> 


»t  -t 


+  I  J  F<rfJk(t-T’t-T')l'ai(T)''.jk(l')dT  dt| 


-00  -00 


„t 


J  J  it  ral( t-T,t-Tl)lCai(T)T(Tl)dT  dTl  , 


(3-5) 


and  we  have  made  use  of  the  symmetry  conditions  (3.5) . 

nie  inequality  (3.*0  must  hold  for  arbitrary  continuous  histories 

•  *  • 

of  y^t),  Kia(T)>  t(t)  •  71:16  integral  operators  in  (3-4)  depend 

smoothly  upon  the  histories,  and  thus  changing  the  histories  in  the 
neighborhood  of  the  present  time  t  =  t  produces  only  a  small  change  in 
the  value  of  the  integrals.  If  the  present  rates  could  be  assigned  in  a 
completely  arbitrary  manner  without  changing  the  values  of  the  coefficients 
of  the  present  rates,  then,  clearly,  in  order  that  the  inequality  be 
satisfied  it  would  be  necessary  that  the  coefficients  vanish  identically. 

It  can  be  shown  that  this  is  also  true  in  the  present  case.  The  formal 
argument  is  similar  to  that  of  Coleman  [  7)  J  for  the  interested  reader 
the  details  are  given  in  an  Appendix. 

The  vanishing  of  the  coefficients  then  yields 


13. 


pt  •  ft 

pS  =  -  x(0)  -  J  ni(0,t-T )t(t)<1t  -  J  $ij(t-T,0)Yij('r)dT 


ft 

'  J  rai(t-T,°)Kai(T)dT  , 


(3.6) 


<*P 


"  2^(0) 


J1 


Gt>eki(0>t-TKi(T)dT  ♦ 


r 


$ae(°,t-T)T(T)dT 


+ 


Yiap 


(t-T,0)Kvl(T)dT] 


(3.7) 


"k 


-  2[Dk3(0) 


J  ak3mn(0>t-T)V(T)dT  +  J 


$k3(0,t-T)T(T)dT 


Faik3(t-T’°)Kai<T)tiT 


5 


(3.8) 


aa 


=  H  .(0) 
ai 


r 


HolBj<0’t'T)l<BJ(T)dT  + 


f 


Foijk(0>t-T)Yjk(T),iT 


-00 


(3.9) 


The  inequality  (3.4)  then  reduces  to 


-J  at  Dij(t*T)'fij(T)dT  ‘I  I  Mw)i(,)dT 
—00  —00 

-  f  lr  H«i(t-T),iai(T)dT  ■ A  -  r  i  s  0 


(3.10) 


14. 


The  first  three  terms  in  (3.10)  are  of  first  order  in  the  rate  histories, 
while  the  remaining  terms  are  of  second  order.  Hence,  to  satisfy  (3.10), 


we  must  have 


,  i3T,n 

T0  *x 


and 


"  J  at  Dij(t-T)YiJ(T)dT  -  J  ^  X(t-T)T(r)dT 

•ffi  —CO 

'  J  &  Hoi(t"r)l'ol(T)dT  s  0  • 


•  •  • 

As  Y..(t)  ,  T(t)  and  k  .  (t)  can  be  chosen  arbitrarily,  we  must 
J  Qfl 

have 


•  0  -  * 0  •  |rHai(t>  * 0  • 

Since  (3 .11)  must  hold  for  homogeneous  temperature  distributions,  we 
also  have 

-  A  S  0 

From  (3.11) »  we  see  that  the  familiar  condition 

-  h  §£  ?  0 

ax 

does  not  necessarily  follow  from  the  entropy  production  inequality  for 
a  viscoelastic  rod.  However,  to  satisfy  (3.1l)  it  is  sufficient  that 
(3.15)  be  valid.  If  we  restrict  our  attention  to  a  class  of  materials 
for  which  (3 .15)  holds,  it  is  easily  seen  thaip  the  constitutive  equation 
(3.3)  must  reduce  to 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


15. 


where  K  (>  0)  is  a  constant.  The  arguments  leading  to  (3*13)  and 
(3.16)  are  similar  to  those  of  Christensen  and  Naghdi  [4]  in  their  work 
on  the  general  linear  viscoelastic  solid. 

Substituting  from  (3*1)  and  (3*6-9)  into  the  energy  equation  (2.17) 
and  neglecting  second  order  terms,  we  obtain  a  linearized  energy  equation 

-  p  To  S  +  pr  -  ||  -  0  (3.17) 

which  can  be  used  to  determine  the  temperature . 


4 .  Symmetries 


We  restrict  our  attention  now  to  a  rod  whose  Helmholtz  function  is 
invariant  under  the  transformations: 


X  -  t  X  ,  A.  -  t  A.  ,  A,  -  i  A, 


where  we  may  take  any  combination  of  +  and  -  .  It  is  a  straightforward, 
but  tedious,  calculation  to  find  the  form  of  the  kernel  functions  under 
these  restrictions .  We  present  here  only  the  final  results .  The  free 
energy  A  is  now  given  by 

Pqa  *  P0A0  "  P0SoT  +  J  J  ^  2  6l(*^12^T^12^  +  2  g2^  ’  ^23^Y23^ 

— OD  —00 

^  83  (  J  *  ^11^)^11^)  ^  ^  S^(  >  )Y22^t )*Vg2 CT( ) 

+  £  86(j)Y33(t)y33(T1)+ ^  87(,)[y11(t)y22('T1)+ Y22(t)y1i(T1)] 

+  8g(  j  )f  Y22 (t  ^33(^1 )  +  Y33 (t )y22(^1)  ^  +  8^( » )(  Yjl^^ ) Y 33 ( *f) )  + 

••  1*  1*  1*" 
YssC'OYiit'n)]  +  ^2  cP1(>)y11(t)+  2  92(»)y22(t)+  2  «P3(  ’  )v33(t)  ]T(Tl) 

+  |  m(,)T(T)T(Tl)+|  h1(,)K11(T)K11(n)+ |  h2(,)[lC11(T)K22(^) 

+  ^22(t )^h (T| ) +  2  ^3 (  s  )^22(*^ )^22(*^)  ^  ^4 (  *  )*1P(T|) 

+  |  h5(  ,)[ki2(t)K21(T|)  +  K21(t)k12(T1)]+  \  h6(,)K21(T)K21Cn) 

+  \  h7(,)K13(T)lC13(Tl)  +  |  hg(,)K23(T)K23(Tl)  }  dT  dT^  , 

where  we  have  used  the  notation 


(4.1) 


(4.2) 
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U 


f(,)  H  f(t-T,  t-Tl)  , 

and  we  have  set  -  \(0)  =  dSq  and  assumed  that  D^O)  =  0  ,  i.e.,  the 
rod  is  free  of  initial  stresses. 

The  constitutive  equations  (3-6),  (3*7),  (3.8)  and  (3*9)  become 


nl  =  J  g3(0,  t-T)v13(T)  dT  , 


=  J  g2(°>  t-T)v23(T)  dT  , 


(4.3) 


= 


AW 

J  [g6(0,t-T)v33(T)  +  gg(05t-T)v22(T)  +  g9(0,t-T)y11(T)]  dT 

—OL' 

-t 

+  J  ^3 (0 5 t“T ) T(t )  dT  , 


pt 

ml  =  J  h8(°5t’T^K23^T^  dT  » 

-00  J 

pt 

m2  =  -  J  hy(0,t-T )k13(t )  dT  ,  (4.4) 

-05 

pt 

ra3  ■  J  (  (hi+-h5)(0,t-T)K12(T)  -  (h6-h5)(0,t-T)K21(T)3  dT  , 

“00 


P11  =  J  ^hi(°’t"T)Kn(T)+  h2(0,t-T)K22(T)]  dT 


A  V  *  * 

P22  ~  J  ( ^3(® ) K22  (t ) +  hg ( 0 , t _T )k^^ ( T ) ]  dT  > 


pt 

P12  +  P21  =  J  [(hU+h5)(°>t-T)^i2(T)  +  (h6+h5)(°5t-T),C2i(T)]  dT 


(4.5) 
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o 


nil  =  J  [S4(°»t-T)V11(T)+  g7(0,t-T)v22(T)+  g9(0,t-T)Y33(T)]  dT 
|»t 

+  J  cp^Ojt-Oi^T)  dT  , 

-00 

pt 

n22  =  J  tg^(0,t-T)v22(T)  g7(0,t-T)vi;L(T)  +  g3(0,t-T)v33(T)]  dT  (4.6) 
ft 

+  J  tp^(0»t-T )T(t)  dT  , 

-00 

rt 

n12  "  n21  =  J  SiCO^-tJv^Ct)  dT 

-00 

and 

pS'  *  pSo  -  m(0,t-T)T(T)dT-|J  [cp1(t-T,0)Yil(T)+ cp2(t-T,0)v22(T) 


+  cp3(t-T,0)Y33(T)]  dT 


(4.7) 


The  quantity  A  is  now  given  by 
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.t 


A  =  J  J  {  Ijt  8xC»)y12(t)y120i) 


(4.8) 


+  2  at  »  )v23(t )v23(ir))  +  2  at  g3^  ’  )v13(T)v22C'n) 

+  1+  at  g2(>)Y1i('r)v1i(Tl)  +  1+  at  g5(»^22(T)V22(Tl) 

+  £  ^  g6(’)Y33(T)^33(ll)  +  k  zt  87(»)v11(t)y2201) 

+  2  at  g8^  ,^y22^t^y33^^  +  2  at  g9^  *^11^^33^ 

+  ^2  at  <pi^»^1i(T)  +  2  at  ^(»^22(t)  +  2  at  t(ti) 

+  2  at  m(  »)t(t)t(ti)  +  2  at  hi(  »)Kii(’r)^ix(T1) 


at  h2^Kll^K22^  +  2  at  h3^  ’^22^*22^ 


I  It  \(.)ii12(T)i12(il)  +  ^ 


+  2  at  h6^  *21^*21^  2  at  **7^ 


+  2  at  h8^»^23^T^23^^  3  dT  d*n 


Inspection  of  the  equations  of  motion  (2.9),  (2.14)  and  (2.15),  the  energy- 
equation  (2.18)  and  the  constitutive  equations  (4.3-7)  shows  that  they  can 
be  separated  into  four  distinct  groups:  two  governing  flexure,  one  governing 
torsion  and  one  governing  longitudinal  motions.  Thermal  effects  and  the 
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energy  equation  need  only  be  considered  with  the  equations  governing 
longitudinal  motions. 


21 


5 •  The  Flexure  Equations 


The  two  sets  of  equations  governing  the  flexure  of  the  rod  are 


3n1  a2ux 

+  pfl * p  ’ 

dm2 

—  ^  nl  -  pi13  *  p0l  —p  -  0 

OXt 


n. 


Pt 

J  g3(0,t-T)v13(T)dT  ,  m2  -  -  J  h7(0,t-T)K13(T)dT 


(5-1) 


and 


sbi3  9U2 

yi3  “  bi3  +  b3i  ’  Ki3  ax  *  b3i  =  aT~  * 


ar^ 

ax 


_2 

a  u2 


iv  +  Pf2  "  p  2 

at 


am1  a  c>23 

W  “  n2  +  pA23  “  pa2  2° 

O  w 


J  g2(°>t-T)v23(T)dT  ,  m1  =  J  hQ(0,t-T)K23(T)dT 


(5.2) 


ab23 

y23  =  b23  +  b32  ’  K23  *  “ax  5  b32  =  3x" 


As  a  simple  application  of  these  equations  we  may  discuss  the  propa¬ 
gation  of  flexural  waves  along  an  infinite  rodt  To  simplify  our  discussion 

+ Throughout  this  Section,  we  restrict  our  attentions  to  the  case  when 
cr2>  p,  g^>  g^  etc.,  are  independent  of  z  . 
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we  introduce  the  notation 


G(s)  =  gg^s)  ,  H(s)  =  hg(0,s) 

Assuming  that  the  waves  are  of  the  forms 
Ug  «  Ue-*^  cos  o)(t  -  , 

b2^  =  B1  e"^x  cos  u)(t  -  ~)  +  B2  e”^x  sin  cu(t  -  x)  , 
on  substitution  of  (5-4)  into  (5. 2)^.  g  we  find 

ng  *  [-0)(B1-hU)Gc  +  0)(B2  +  |  U)Gs]e"^X  sin  u>(t-^) 

+  [«)(B1-jjiU)Gs+ tu(B2+^  U)Gc]e‘^X  cos  m(t  -  ^ )  , 

m x  -  [u)(nB1+“  B2)Hc  +  u)(-|iB2  +  “  B-^H^e  M'X  sin  co(t-^) 


(5.3) 


(5.4) 


(5.5) 


+  [ “U)(m<B1  +  “  B2)Hc  +  U)(-nB2  +  |  B1)Hs]e"^X  cos  u>(t  -  £)  , 

where 

GB(«)  -  r  G(u)  sin  o)U  du  ,  Gc(u>)  =  J  G(u)  cos  cuu  du  ,  (5.6) 

o  o 

and  H  (to )  and  H  (to)  are  similarly  defined  as  the  Fourier  sine  and 
cosine  transforms  of  H(u)  .  Using  (5-5) >  (5-4)  and  the  equations  of 
motion  (5*2^  g>  we  obtain  a  set  of  homogeneous  equations  to  determine 
U,  B1,  and  Bg 
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t-2|*  f  «,+  (n.2  -  25)GS  t  CluJu  +  [-»Gs  +  f  GC]B1  -  tnGc  +  f  Gs1b2  =  0  , 


[211|as+((12-%)G=]u-  K*'st]J1*t-*ocnc!)B2.  0  , 


[^s  -  ?  °clu  *  [^Hs  -  2“  !  Hc  *  %  Hs  -  °S  +  P"2“]B1 

c 


*  [2|*fVu2Hc-%  hc-gc]b2  -  0  , 

c 


[g,G  -  —  G  ]U  +  [-2(j,  —  H  +  —  H  -  \x2  H  -  G  )b, 
^  s  c  c  Kcsc  C  ^  c  sl 


(5.7) 


+  [-  \  Hs  +  n2  H  -2m,  Hc-  Gc+  pa2a)]B2  =  0  . 

c 


These  equations  will  determine  ,  c  and  the  direction  of  the 
vector  {U  ,  B.^  ,  B2)  as  functions  of  the  frequency  cu  .  The  solution 
for  an  arbitrary  value  of  the  frequency  is  quite  complex,  and  we  will 
restrict  our  attention  to  the  asymptotic  behavior  of  the  solution  as  the 
frequency  becomes  very  large.  If  G(u)  and  H(u)  are  sufficiently 
smooth,  it  is  well  known  (see,  e.g.,  Lighthill  [8])  that 


Gs^ 


0)  3 

II) 


Gc(tu)  ~ 


•  •  • 


(5.8) 


with  similar  expansions  for  Hs(u))  and  Hc(u))  .  We  assume  that  b  , 
Bj_  »  Bg  ,  p,  ,  and  c  have  asymptotic  expansions  of  the  form 
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> 


> 


U  ~  u 

o 
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B5  ~  b9n  +  -Si  + 

2  20  U) 


>*1 


U,  ~  LL  +  -  + 

P  pO  U) 


(5.9) 


c  ~  c  +  — 
O  (I) 


Substituting  (5.8)  and  (5-9)  into  (5-7)j  and  then  investigating  the  con 
ditions  that  non-trivial  solutions  exist,  we  find  two  modes  of  propaga¬ 
tion  for  flexural  waves .  In  the  first , 


c 


2 

o 


G M 

P 


uo  4  0 


10 


20 


(5.10) 


while  in  the  second 


c2  »  HM 

o  p  a2 


^o 


> 


u 


10 


20 


4  o 


(5.11) 


These  relations  generalize  well  known  results  for  extensional  waves  in  a 
one  dimensional  visco-elastic  medium  (see,  e.g.,  Hunter  [5  ))• 

Our  developments  here  rest  on  the  initial  postulates  of  a  balance 
of  energy  and  an  entropy  inequality.  For  a  curve  with  assigned  directors 
embedded  in  a  Euclidean  3-space  the  resulting  equations  in  Secs.  2  and 
3  are  exact  and  aside  from  linearization  involve  no  approximations. 
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However,  the  equations  given  here  could  be  obtained  from  a  three  dimen¬ 
sional  theory  of  rods  after  suitable  approximations  are  made.  A  discus¬ 
sion  of  the  kind  of  approximation  necessary  is  given  in  [3]. 

To  one  familiar  with  classical  theories  of  rods,  the  form  of  our 
equations  is  quite  suggestive,  and  we  will  present  here  a  few  rather 
intuitive  ideas  concerning  the  connection  between  the  two  theories. 

A  rod  as  a  three  dimensional  body  may  be  regarded  as  the  Cartesian 
product  of  its  line  of  centroids  and  its  cross  section.  Because  a  rod 
is  a  "thin"  body,  we  can  constrain  the  cross  section  to  undergo  only 
homo^neous  deformations .  The  homogeneous  deformation  of  the  cross 
section  is  then  determined  by  the  motion  of  two  independent  vectors. 

Thus  in  the  model  of  the  Cosserat  curve,  the  two  directors  can  serve  to 
describe  the  homogeneous  deformations  of  the  cross  section.  Since  any 
two  independent  vectors  uniquely  determine  a  homogeneous  deformation,  we 
may  choose  to  take  the  directors  along  the  principal  axes  of  inertia  of 
the  cross  section. 

These  ideas  will  be  sufficient  to  permit  us  to  compare  the  results  of 
the  Cosserat  theory  with  solutions  from  the  theory  of  linear  viscoelasticity 
for  an  isotropic  three  dimensional  rod.  We  consider  first  the  case  of  iso¬ 
thermal,  quasi-static  pure  flexure.  Using  the  usual  quasi-static  "corre¬ 
spondence  principle"  between  elastic  and  viscoelastic  solutions,  the  solution 
of  the  viscoelastic  problem  in  the  Laplace  transform  plane  is  easily  obtained 
from  the  corresponding  elastic  solution.  We  consider  a  straight,  isotropic 
prismatic  rod  and  introduce  a  Cartesian  coordinate  system  (x  ,  y  ,  z) 
where  the  z  axis  coincides  with  the  line  of  centroids  of  the  cross  section, 
and  the  x  and  y  axes  are  taken  along  the  principal  axes  of  the  cross  section. 
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Consider  now  the  problem  of  pure  flexure  of  such  a  body.  If  a  denotes  the 
transform  of  the  stress  tensor,  the  only  non-zero  component  of  the  stress  is 


CT33  =  ax  +  By  , 


(5-12) 


where  a  and  3  are  constants.  The  three  dimensional  displacement 

* 

vector  is  denoted  by  u  and  its  transform  is  given  by 


ux  =  -  ~  (|ax2  +  pxy  -  ±  &yc)  - 


1  „.2v  a  2 


sE 


2sE 


-*  v  .  12  1  0  2.  B 

Up  =  -  —  (Qfxy  +  p  0y  ■  j  gx  )  - 

sE  2sE 


(5.13) 


Up  =  *“  (<*x  +  0y)  z  , 
0  sE 


where  v  and  E  are  defined  in  terms  of  G^t)  ,  the  relaxation  func¬ 
tion  in  shear,  and  Gg(t)  ,  the  relaxation  function  in  isotropic  com¬ 
pression,  by 


X  2g2  +  Gi 

E 


v  = 


G2  '  G1 
2G2  +  G1 


(5.14) 


The  transforms  of  the  moments  acting  across  any  section  are  then 


m,  =*  31 


xx 


sEI 


-2—* 
d  ug 

xx  .  2 
Sz 


x,y=0 


(5-15) 


m_  -  -  Ofl  *  sEI 

2  yy 


J2-* 
5  U1 


x,y*0 
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where  I  and  I  are  moments  of  inertia  about  the  x  and  y 
xx  yy  J 

respectively. 

Taking  the  transforms  of  (5.1)  and  (5.2)+  we  easily  obtain  the 
equations 


axes 


n1  =  0  ,  nig  =  -  sh^ 


3b 


±1 


az 


and 


-  -  -  au: 
b13  +  b31  '  0  >  b31  *  5T 


_  ab 


=  0  ,  mx  -  shg  — 
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au2 

b23  *  b32  ’  0  ’  b32  *  5T 


Thus 


-2— 

_  _  a  u2 

=  -  shg  — , 

°  az 


-2— 
a  u. 


m?  =  Sft7  — 2“  » 

f  az 


(5.16) 


(5.17) 


(5.18) 


and  comparing  these  results  with  (5.15)>  we  obtain 


8 


El. 


xx 


(5.19) 


We  assume 


^(t) 


and 


bij(T) 


vanish  on  -  »  <  t 


0  . 
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6.  Torsion 


The  equations  governing  torsional  motion  of  the  rod  are 


dnu  d2t>P1 

&z  +  p^12 "  l21^  “  p(al  2  “  a2  2  ^  * 

ot  ot 


a(p12  *  P21) 

az 


2  2 
3  b12  *  b21 

+  pU12+  l21)  ~  2  +  °2  2  ^  =  2tt12 


pt 

m3*J  ^h4"h5^0,t“T^Ki2^T^  “  (h6'h5^°»t‘T^21^T^dT  » 

—08 

.t 

P12  +  P21*J  t(h4+h5)(°»t“T)^i2^T^+ (h6+h5^0,t‘T^21^T^dT  » 


(6.1) 


rb 

n12  *  J  «i(0»t“T)Y12(T)dT  , 


^b12  ^b21 

K12  *  —  *  K21  “  —  *  V12  =  b12  +  b21 


We  restrict  our  attention  here  to  the  case  when  the  free  energy  and 
kinetic  energy  are  invariant  under  the  transformations 


£l  “*  £2  ’  £2  "*  £ 1 


(6.2) 


This  restriction  implies  that 


«1  “  ®2  »  h4  *  h6  * 


(6.3) 


Under  these  conditions,  we  seek  solutions  of  (6.1)  with 
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0 


> 


(6.4) 


'12 


12 


bj>) 


Then,  if  the  assigned  loads  are  assumed  to  vanish,  equations  (6.1) 
reduce  to 


where 


^3  2  *  iu>t 

~  P«“  bi2  e 


_  ?  6  db12  iu)t 

m3  =  2  H1  ~  e 


0(  ■  2<y^  =  2q»2  , 


^(Oju)  =  h4(0,u)  -  h5(0,u)  , 

H  =  itu  f"  H.(0,u)  e"iuni  du  , 
1  J0 

-  its  ^(Ojico) 


Combining  (6.5  a,b),  we  obtain 


2  * 

d  b12  2 


where 


T  • 


(6.5) 


(6.6) 


(6.7) 


(6.8) 
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If  we  take  b^g(O)  =  0  ,  the  solution  of  (6.7)  is 


b12  =  b  sin  q  z 


(6.9) 


The  torque  necessary  to  maintain  this  deformation  is 


m„  =  b 


POU) 


cos  q  z  e 


iu>t 


(6.10) 


This  result  may  be  compared  with  Berry's  [  6  ]  result  for  the  torsional 
vibrations  of  a  circular  cylinder  obtained  from  the  three  dimensional 
theory.  With  some  changes  in  notation.  Berry ' s  expression  for  the  torque 
necessary  to  maintain  sinusoidal  oscillation  is 


2  2 

m3  =  b  P2X^  C°S  X  Z 


x2  = 


2pti>£ 


■na2[iu)G1(iu))] 


(6.11) 


where  a  is  the  radius  of  the  cylinder.  The  tangential  displacement  of 
points  in  the  cross  section  is  given  by 


u  =  b  r  sin  q  z  e 
0 


ittJt 


(6.12) 


Thus,  b^  can  be  interpreted  as  the  rotation  of  a  radial  line  element 
in  the  cross  section.  Equations  (6.11)  and  (6.10)  are  equivalent  if 


a  -  %■  >  ^(i®)  =  S  Jz  Viu>) 


(6.13) 


where  Iz  is  the  polar  moment  of  inertia  of  the  cross  section. 
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7.  Longitudinal  Motions 


In  the  theory  employed  in  this  paper,  thermal  effects  only  arise  in 
the  case  of  extensional  motion  of  the  rod.  We  restrict  our  attention  to 
the  case  when  the  heat  flux  is  given  by  (3.19)*  Then,  the  equations 
governing  the  motion  are 


^3 

_  j-* 

2 

a  u 

ax 

+  px, » p 

y 

nll 

3pu 
ax  + 

p(iu  - 

“i 

a\i 

at2 

)  , 

OJ 

OJ 

Sp22 
ax  + 

P ( ^22  “ 

*2 

at2 

)  , 

-  p 

Tq  S  +  pr 

ah 
“  ax 

0 

y 

n3  = 

pt 

J 

-CD 

-t)v33(' 

r)  4 

g8(°. 

•  • 
t-T)Y22(T)  +  g9(°>t“T)Y11(T) 

+  <p3(0,t-T)T(T)]dT  , 

l*t 

pH=J  [h1(0,t-«r)K11(T)  +  h2(0,t-T)K22(T)]dT  , 

-CD  '  ' 

»t 

p22  =  J  [h3(0>t-T)K22('r)  +  h2(0,t-T)K11('r)]dT  , 

-CD  ^ 

TTH  =  J  [g4(0,t-T);il(T)  +  gyCOjt-T^gtT)  +  g^Ojt-lOv^T) 


(7.1) 


+  ^(O.t-OTC'OW  , 


32. 


,t 

n22  “  J  +  S7(°»t-T)V11(T)  +  g8(0,t-T)v33(T) 


+  q^(0,t-T)T(<r)]dT  , 


Yu  -  2b^  ,  y00  m  2b, 


du. 

-  2  r— * 


22  “‘*‘'22  *  Y33"Z5t“  *  K11 


!^u  „ 

3x  *  *22  “  fcx 
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Appendix 


The  inequality  (3.4)  can  be  written  in  a  more  revealing  form  if  we 
introduce  the  vectors 


3a  (2  Yy.  T) 


(Al) 


£  -  ("ofl’  V  *0fi»  ■  pS/) 


Then  the  constitutive  hypothesis  for  n.  ,  p  rr  „  and  S#  in  terms 

X  QfX  CmP 

of  a  matrix  kernel  function  P(t)  and  a  constant  vector  T  is 


£  •  Eo  +  J  P(t-T>3(T)dT  . 


(A2) 


By  introducing  appropriate  matrix  kernal  functions  g(s,t)  and 
H(t)  in  (3-1)  we  can  write  pA  in  the  form: 


•  pi 

pA  »  [h(0)  +  G(0,t-T)^(T)dT]  ^(t) 

~  ~  ~  ^ 


t  „t 


J  »  fi(t‘T)2(T)dT  *  J  I  £  £(t-T,t-s) 


l.t  «t 


J  I  ^  GCt-T^-s^TMsJdTds 


(A3) 


The  inequality  (3*4)  for  processes  with  ~  (t)  =  0  then  becomes 

Oz 
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rt 

(r  -  H(0)  +  [p(t-T)  -  G(0,t-T)]ll(T)dT}  n(t) 

~  ~  «  A*  ~  M  M 


»t  n't 


-J  ^  H(t-T)Tj(T)dT-  J  J  ^  G(t-T,t-s)Tj(T)T](s)dTds  S  o  .  (A4) 


We  will  assume  that  T1(t),  ~  H(t),  -rr  G(t,s),  and  —  G(t,s)  are  con- 

^  Ob  r*  Ot  **  OS  ^ 

tinuous  (these  assumptions  have  already  been  used  when  differentiating 
under  the  integral  sign  to  find  pA)  .  Typical  terms  in  (a4)  are  of  the 
form 


J  Cp(t-T)'n(t)dT  , 

i|i(t-T,t-s)Tl(T)T](s)dTds  , 


ft 


(A5) 


where  cp  and  i(r  are  continuous  functions.  If  ^(t)  is  some  specified, 

continuous  history  and  if  a  is  an  arbitrary  vector,  then  given  any 

A.  .  A 

e  >  0  we  can  find  a  continuous  history  T|(t)  such  that  T|(t;  =  a  and 


||  J  tp(t-T)Tl(T)dT  -  J  ^(t-T)ll(T)dT||  <  e  , 

V  n/  v  A#  rss 

-00  -00 


II  r  f  ♦(t-T,t-s)‘n(T),n(s)dTdS  -  J  J  t(t-T,t-s)'n(T)Tl(s)dTds|| 


<  e 


The  norm  in  (A6)  is  the  usual  vector  or  tensor  norm;  i.e.,  if  V  is  a 
vector  l|v||  ■  [v  v]1^2,  if  A  is  a  tensor,  then  ||a||  =  tr[A  a^]^2 . 

■W 1  /sw  /W  r*j 

For  example,  consider 


(A6) 


ti(t)  =  TKt) 


-eo  <  |  S  t  ■  6 


=  Tl( t—  6 )  +  [of- Tj(t-6)3  [— .  t  -  6  S  t  S  t  (A7) 
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Then 


||  J  cp(t-T)T|(T)dT  -  J  cp(t— t )T)(*t )d*T"J| 

*  Arf  A#  W  ^  /W 

—00  —00 

=  ||  J  cp(t-T)£*n(T)  -  T](t-6)-  [a- Tl(t-6)[I~^]}dT||s$(3r  +  ||a||)6  ,  (A8) 

t-6 


where  $  =  sup||<p(t— t ) ||  and  T=  sup||t\(t)||  on  some  sufficiently  large 

/•w  f*W 

interval  containing  (t-6,t)  .  Similarly 


-00  -00 


*  y1t2  +  (2r  +  |H|)2]62 


(A9) 


where  Y  =  sup||^(t-T,t-s)||  .  By  choosing  6  sufficiently  small  we 
obtain  our  "desired  result. 

Consider  now  the  inequality  (A4)  for  a  given  history  T|(t)  .  Choose 

A  A  ~ 

another  history  T1(t)  with  T|(t )  =  a  ,  a  being  an  arbitrary  vector. 

<v  ~  IS/ 

This  history  must  also  satisfy  the  inequality  so  that 

(•"t  A 

(r  -  H(o)  +  [ P(t~T )  -  G(0,t-T)]*n(T)dT}  a 

-00 

-J  ^  H(t-x)j(T)dT-J  J  ^  G(t-T,t-s)T}(T)T](s)dTds  ?  0  .  (A10) 

—00  “00  —00  E 

Adding  and  subtracting  terms  in  (A10)  we  obtain: 
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p't  A 

{  I  C  P(t“T )  -  G(0,t-T)]  C^(t)  -  *n(T)]dT}  •  a 

*  „  ^  /V  ^ 

“00 

-J  ^  H(t-T)[j(T)-Tl(T)]dT 


»t  »t 


n^r  G(t-T,t-s)C^(T)Tl(s)  -  Tl(T)Tl(s)]dTds 

•03  “00 

{r  -  H(0)  +  [p(t-T)  -  G(0,t-T)]‘n(T)dTl  •  a 

s-*jO  rv  J  /-w  r*t 


J  ^  H(t-T)Tj(T)dT  -  ^  f  ~  G(t-T,t-s)Tj(T)Tj(s)dTds  *  0  .  (All) 


however,  using  (a6)  we  see  that  by  choosing  Tl( T )  properly  the  first 
three  terms  in  (All)  can  be  made  arbitrarily  small,  and  hence  if  the 
inequality  is  to  hold  we  must  have 


-00 


+  J  ^  H(t-T)lfj(T)dT  +  J  J  ^  G(t-T,t-s)Tj(T)Tl(s)dTds  *  0 


(A12) 


Since  a  is  arbitrary  and  T1(t)  is  any  continuous  history,  (A12)  implies 
that 


(A13) 


and 


P(t-T)  =  G(0,t-T) 


(Al4) 
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The  inequality  (a4)  then  reduces  to 


-J  4  H(t-T)T](T)dT  -  J  J  ^  G(t-T,t-s)Tl(T)Tl(s)dTds  5  0  .  (A15) 

-00  ~  -00  -00  °  ~  ~ 

* 

Let  T]  (t)=  P  T1(t)  where  p  is  an  arbitrary  scalar.  The  inequality 
(A15)  must  still  be  satisifed  for  this  new  history  so  that 


-p|  |  H(t-T)T)(T)dT  -  P2  jrt  f  £  G(t-T,t-s)Tj(T)3(s)dTds  £  0  .  (Al6) 


In  order  that  (Al6)  be  satisfied,  we  must  have 


^  H(t-T )  =  0  ,  (A17) 

and 

G(t-T  ,t-s)T|(T)T|(s)dTds  >  0  .  (Al8) 

ou  ~  ^ 

Our  original  inequality  (3.4)  in  view  of  (A13),  (Al4),  and  (A17)  reduces 
to 


rr 
-00  -00 


_h_  ST 
Tq  Sz  " 


0 


(A19) 


If  we  revert  to  our  usual  index  notation,  (A13)  and  (Al4)  imply  (3*6), 
(3-7),  (3.8),  (3.9),  and  (3.10). 
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